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Interacting classical and quantum particles 
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We apply Hall and Reginatto's theory of interacting classical and quantum ensembles to harmon- 
ically coupled particles, with a view to understanding its experimental implications. This hybrid 
theory has no free parameters and makes distinctive predictions that should allow it to be ex- 
perimentally distinguished from quantum mechanics. It also bears on the questions of quantum 
measurement and quantum gravity. 
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Does physics contain a classical sector, such as gravity? 
How does the quantum sector react back on the measure- 
ment devices that collapse it? What is the correct theo- 
retical description of nano-mechanical systems in which 
quantum and classical objects are coupled? Such ques- 
tions motivate the investigation of theories of interacting 
quantum and classical systems. 

A number of candidate theories have been developed, 
but most have been found wanting in some aspect [J- 
3]. In this paper we investigate the ensemble theory of 
Hall and Reginatto 041], which avoids the difficulties in 
certain other attempts to unite quantum and classical 
dynamics within a unified mathematical structure. 

Experimental investigation of objects straddling the 
boundary between the quantum and the classical is be- 
coming possible with nano-mechanical techniques 043, 
such as coupling of cantilevers or membranes to Bosc- 
Einstein condensates [l(J, and with electromagnetic sys- 
tems, such as SQUID rings coupled to classical circuits 
[ill ]. In so far as this boundary is experimentally un- 
explored, a quantum mechanical description of such ex- 
periments [12J should be compared against the foil of an 
alternative theory. 

Amongst the best developed alternatives to quantum 
mechanics are the collapse theories that postulate a tran- 
sition between quantum and classical mechanics at a scale 
determined by a free parameter [H, I3- In contrast, the 
hybrid theory of Hall and Reginatto postulates distinct 
quantum and classical sectors, not necessarily a transi- 
tion between them. It is a different kind of theory to 
collapse theories, and in particular it has no free param- 
eters. Its predictions may be compared with those of 
quantum mechanics and of collapse theories. 

There are difficulties with a standard quantum field 
theory of gravity [l5j , and a number of lines of research 
are seeking a consistent quantisation method. Although 
there are good reasons to want to quantise grav ity |16j . 
classical gravity has not yet been ruled out [17h20| . de- 
spite one par ticular theory having been experimentally 
disproved [21j . One of the obstacles to embracing classi- 
cal gravity is the lack of well-developed theories encom- 
passing both the quantum and classical sectors. The hy- 
brid theory of Hall and Reginatto achieves this and has 
been applied to the case of a classical gravitational wave 



interacting with a quantized scalar field [22| , providing a 
counterexample to claims that this is not possible [231 . 124] . 

Some formulations of quantum mechanics require a 
classical sector to accommodate measurement devices. 
The lack of a quantitative theoretical distinction between 
the classical and quantum sectors is one aspect of the 
quantum measurement problem [25]. Indeed, axiomati- 
sations of quantum mechanics usually include measure- 
ments as an unanalyzed notion [26]. Another aspect of 
the measurement problem is that within standard theory, 
classical systems act on quantum systems, by reduction 
of the state, but quantum systems do not react back on 
classical systems. This back action is normally negligible 
when the measurement apparatus is large compared to 
the quantum system. However, at the quantum-classical 
boundary this may no longer be true. 

Besides their fundamental interest, theories of inter- 
acting classical and quantum systems may be useful for 
developing consistent approximations of subsystems as 
classical, for example in quantum chemistry [271 ] . 

In the following pages, we apply the theory of Hall and 
Reginatto, which we shall refer to as the hybrid theory, 
to a classical particle interacting via a harmonic poten- 
tial with a quantum particle, Fig. [1] Our purpose is to 
understand the predictions of the theory and how they 
compare to the familiar cases of two interacting classical 
particles or two interacting quantum particles. In partic- 
ular, we investigate the coupling between the relative and 
centre of mass degrees of freedom that is characteristic of 
the hybrid theory 0, Q , and find distinctive experimental 
signatures of hybrid dynamics. 

We now outline the hybrid theory. A full account may 
be found in the papers of Hall and Reginatto H-OJ- In 
this paper, we restrict our attention to the simplest case 
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FIG. 1: (color online) Schematic illustration of the physical 
system considered in this paper. The spheres represent the 
particles, and the arrow the single dimension to which they 
are confined. The dashed line between the particles represents 
their harmonic interaction. Otherwise, they are free to move. 
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of one spatial dimension, with one classical particle and 
one quantum particle. Let the quantum coordinate be 
q and the classical coordinate x. The theory is based 
on probability densities in configuration space. Let the 
probability density for finding the quantum particle at 
q and the classical particle at x be P(q,x). Specifying 
the dynamics of P{q, x) in the Hamiltonian formalism 
requires another function S(q, x) that is canonically con- 
jugate to P{q, x), and an ensemble Hamiltonian 



H[P,S] = dqdxP 



(d q S) 2 , (d x S) 2 



2m„ 
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+V(q,x,t) + 
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The masses of the quantum and classical particles are m q 
and m x respectively. The interaction potential between 
the particles is V(q,x,t). We define the Poisson bracket 
[El. [2a] between two functionals A[P, S) and B[P, S] as 



{A,B} 



dq dx 



8 A SB 
~5P~5S 
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For a functional L[f] = J dzF(f, d z f, z), for example, the 
variational derivative has the usual form [281 ] 



SL/8f = dF/df-d z [dF/d(dJ)]. 



(3) 



Using this, the dynamical equation for any functional 
A[P, S, t] is i 

dA/dt = {A, H} + dA/dt. (4) 
In particular, the dynamics of P and S are given by 



dP 
~dt 
dS_ 
~dt 
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2m. 



d q {Pd q S) 
■{d q Sf 



-d x (P8 X S) , (5) 
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■{d x S) 2 -V 
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~2m q P 1 / 2 ' 
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The final term may be recognized as the quantum poten- 
tial [29j. These equations may be combined by defining 
the complex function ip(q, x) = y/P(q, x) e lS ^ q ' x ^ >i . Eqs. 
(I5I6|) are then equivalent to 
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~dt 
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The first three terms on the right hand side of this equa- 
tion have the familiar Schrodinger equation form. The 
final, nonlinear, term ensures that the particle with co- 
ordinate x is classical 130-331. 



Observables in the hybrid theory are defined to be 
as similar as possible to those in quantum and classi- 
cal physics. They are ensemble expectation values, de- 
fined by straightforward extensions of the usual defini- 
tions H3]- For any real classical phase space function 
f(x,p x ), the corresponding expectation value is 



dq dx P /(a;, d x S). 



(8) 



For any Hermitian operator M associated with the quan- 
tum particle, the corresponding expectation value is 



(M) 



dq dxip*Mip. 



(9) 



Expectation values of products of classical and quantum 
observables may be defined based on various principles; 
for example, the exact uncertainty principle [35| . We de- 
fine them so that various straightforward principles give 
consistent results. A simple prescription is to first pro- 
mote classical quantities to quantum ones, use the quan- 
tum definition Eq. ©, and then take the classical limit 
of the promoted quantity by setting h — > only where 
it originates from the classical observable [4]. The only 
such products we will use are 



(qx) 



and 



(PqPx 



dq dx P q x, 



dq dx P d q S d x S, 



(10) 



(11) 



which have the same form as if both particles were clas- 
sical. The second also follows from the symmetrized ob- 
servable \{p x p q +p q p x ) by the straightforward procedure 
of using the classical form d x S for p x and the quantum 
operator —ihd q for p q . 

The hybrid theory has previously been applied to find 
analytic expressions for coherent states of the hybrid har- 
monic oscillator with the interaction potential [J] 



V(q,x,t) = -k(q- 



(12) 



In this paper, we extend that analysis to the more gen- 
eral squeezed Gaussian states. Hence their covariances 
become time-dependent rather than static. By mak- 
ing ansatzes, we reduce the partial differential equa- 
tions Eqs. (|5|6[) to ordinary differential equations that are 
straightforward to solve numerically. 

We formulate the problem in matrix- vector form with 
vector coordinates £ = (q,x) T [4j. The dynamical equa- 
tions ([5]) and (O are then 



dP 
~dt 
dS_ 

~dt 



V ■ (PUVS) = 0, 

i(VS) T (7VS + i£ T C£ 
h 2 V ■ {EiUEiVVP) 
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(13) 
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where V = (d q ,d x ) T , and we have used the matrices 



C = k 



1 -1 
-1 1 



, U = 




(15) 



Depending on the choice of the matrix this form of 
the equations applies to either the classical E c , quantum 
E q , or hybrid Eh case: 



E h 



1 



1 
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We make the Gaussian ansatz for P: 








2tt 



-e 2 



(16) 



(17) 



where K (the inverse of the covariance matrix) is a 
positive-definite symmetric matrix, \K\ is its determi- 
nant, and a is a vector. This is consistent with the equa- 
tions of motion if and only if S is at most quadratic in 
the coordinates £: 

S% t) = 1(£ - a) T L(d - a )+(3 T (^a) + a, (18) 

where L is a symmetric matrix, /3 is a vector and a is 
a scalar. Substituting these ansatzes into Eqs. f| 13[) and 
([T4")) yields a system of ordinary differential equations [1] . 
The terms quadratic in £ — a yield [36| : 



if f/£ + Lf/if = 0, 



dL h 2 

— + LUL + C- —KEiUEiK = 0, 

at 4 



(19) 
(20) 



where the equations have been written to emphasise the 
symmetry of the matrices. The terms linear in £ — a 
yield: 

*£-Uf} = 0, ^ + Ca = 0. (21) 

These equations are the same as for a classical oscillator, 
and have the solution 



a 



(a + bt) 

CfJ, COS (u!fj_t + i 

/ rn \ 

P = M J 1 - CM w M sin (w M t + . 



-1 



(22) 
(23) 



where a, b, c and (f> are determined by the initial condi- 
tions. — m q m x /(m q + m x ) is the reduced mass, and 
cjp = y/k/ n is the oscillator frequency. Hence the prob- 
ability distribution P is centred on the classical motion, 
as in the quantum case. The zeroth-order terms in £ — a 
yield a differential equation for a, which plays no further 
role in our analysis. 



The differences between the probability distribution 
dynamics of the classical, the quantum, and the hybrid 
cases are found in the covariance matrix Z = K~ x . The 
elements of Z are related to the coordinate covariances 
by: 



Var(g), Z x 
Cov(q, x) — 



c = Var(x), 

(qx) - (q){%)- 



(24) 



We have numerically solved the nonlinear differential 
equations (|19I20[) for the matrices K and L (37} • For our 
initial conditions, we assume that the two particles are 
uncorrelated, such that K and the covariance matrix Z 
are initially diagonal. The momentum covariances may 
be expressed as 



Cov(pj,p k ) = (pjPk) - (p 3 )(pk) 



= LZL 



-EiKEi 



(25) 



In the quantum case, we make S initially linear in the 
coordinates so that the matrix L = 0. The quantum 
particles are then in a minimum-uncertainty state. In 
the hybrid case, we give the classical particle the same 
momentum variance as the quantum particle by taking: 

WO) = 0-5hK xx , L M (0) = Wo) = °- We do not 

present solutions for the classical case as the solutions 
for K and L are singular at a quarter-period, although 
there are no singularities in the observables. Our present 
approach is not suited to the classical case, although it 
works well for the hybrid and quantum cases. 

Fig. [5] shows the time dependence of the error el- 
lipses associated with the position covariance matrix Z, 
for both hybrid and quantum cases, for equal masses 
m q = m x . These ellipses have their axes in the directions 
of the unit eigenvectors of Z, with the semi-axis lengths 
equal to the square roots of the corresponding eigenval- 
ues. Hence the ellipses' orientations may be interpreted 
as giving the directions of maximum and minimum co- 
ordinate variance, with the semi- axis lengths giving the 
values of the square roots of those variances. 

The covariance ellipses in Fig. [2] have similar orienta- 
tions, with the quantum ellipses strictly aligned with the 
center-of-mass and relative coordinate directions: 



R = 



m q q + m x x 
M ' 



r = q — x, 



(26) 



where M — m q + m x is the total mass. Prominent shared 
features are: variance oscillations in the relative direc- 
tion, and variance growth in the center-of-mass direction. 
However, a significant difference between the two cases 
is the rocking of the orientation of the hybrid ellipses, in 
contrast to the fixed alignment of the quantum ellipses. 
This rocking arises from a coupling between the center- 
of-mass and relative coordinates, which underpins the 
signatures of hybrid dynamics considered in this paper. 

We next consider the solution in center-of-mass R and 
relative r coordinates. We will find that the center-of- 
mass has distinctive dynamics in the hybrid case. The 




FIG. 2: (color online) Error ellipses for the covariance ma- 
trix Z, shown at quarter-period intervals over two oscillation 
periods. The solid black ellipses are for the hybrid case and 
the dashed blue ellipses are for the quantum case. The axes 
are the dimensionless coordinates: q = q/xo and x — x/xo, 
where xq is the harmonic oscillator length xo = (ft//ioj fl ) 1,/2 . 
First row: t = 0,T/4,T/2. Second row: t = 3T/4, T, 5T/4. 
Third row: t = 3T/2, 7T/4, IT. Here T = 2tt/o; m is the os- 
cillator period. The initial conditions are given in the text, 
with K xx — Kgg — 1/xq. In units of the reduced mass fi, the 
masses are: m, = m x — 2/i. 



associated momenta are the usual: 



PR = Pq + Px 



Pr 



M 



Pq 



mq 
M 



Px- 



(27) 



Denning the new coordinate vector £' = (R,r) T , the co- 
ordinate transformation matrix T such that £ = T£' is 



T = 



m x /M 
-m q /M 



(28) 



In the center-of-mass and relative coordinates, the corre- 
lation matrix Z is given by 

Z' = (T T KT)~ X . (29) 

Fig. [3] shows the components of Z 1 for the equal mass 
case of Fig. [2j It shows that the spreading is occurring 
in the centre-of-mass coordinate (solid lines), while the 
relative coordinate variance is oscillating (dashed lines). 
For the quantum case, the centre-of-mass variance has 
the quadratic dependence on time characteristic of a free 
particle [38| . confirming our methodology. The most no- 
table difference between the hybrid and quantum cases 
is the behaviour of the correlation between the relative 
and ccntcr-of-mass coordinates Z' Rr . In the quantum 
case, the coordinates decouple and there is no correla- 
tion: Z' Rr (t) = 0; in the hybrid case they are coupled, 
and Fig. [3] shows a growing oscillation. 



FIG. 3: (color online) Coordinate covariances in center- 
of-mass and relative coordinates versus dimensionless time 
t = tuj^ for hybrid (black curves) and quantum (blue curves) 
oscillators. The solution is the same as in Fig. [21 Solid lines: 
Var(_R) = Z' RR . Dashed lines: Var(r) = Z' rr . Dotted line: 
Cov(_R, r) = Z' Br . In the quantum case, Cov(_R, r) = 0. 



Such coupling has been previously noted, as the ensem- 
ble Hamiltonian Eq. (|TJ) in center-of-mass and relative 
coordinates is |j| 



H[P,S] = J dRdrP 



(d R Sf 



2M 



q (d R \nP) 2 

(d r S) 2 



8M 2 
dRdrP 



AM 



dRdr 



2fj, 

(d R P)(d r P) 
P 



+ —(d r \nP) 2 + V(r) 
8m 



(30) 



In particular, the last term couples the center-of-mass 
and relative coordinates. 

The center-of-mass kinetic energy E R , relative kinetic 
energy E r and the interaction energy (V) are given by: 



Er 

(P? 
2fi 



<P 2 



Rl 



1 



2M 2M 

1 m 
2~Mm. 



^=-^ = T77 — (Pl) + ^ 2 



(p 2 x )+2(p qPx )), (31) 
2{ PqPx )), (32) 



{V) = \{r 2 )= k -{{q-xf). (33) 



The total energy E = Er + E r + (V) and the center-of- 
mass momentum (p R ) are conserved, as for the classical 
and quantum cases. The classical oscillator parts of the 
solution, Eqs. (|22|23|) . contribute \Mb 2 to the center- 
of-mass kinetic energy, |c 2 sin 2 (o; M i + <f)) to the relative 
kinetic energy, and |c 2 cos 2 (w Al t + <p) to the interaction 
energy. The sum of all these classical energies is inde- 
pendently conserved. That our numerical solutions of 
Eqs. (I19I20P conserve all the required observables in both 
quantum and hybrid cases is a check on their validity. 

In both the classical and quantum cases the internal 
energy, Ej = E r + (V), and the center-of-mass kinetic 
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FIG. 4: Hybrid oscillator center-of-mass kinetic energy Er 
(solid line) and internal energy Ei = E r + (V) (dashed line) 
versus dimensionless time t = tui^. Masses are equal, m q — 
m x . The total energy is E — 1.125ficj^.The vertical dotted 
lines are at quarter-periods. The contribution of the classical 
motion to the internal energy has not been included. The 
solution is the same as in Fig. [2] In the case of quantum or 
classical physics En and Ej are constants. 



energy, Er, are also separately conserved. In the hybrid 
case this is no longer true, as can be seen in Fig. flU). In 
this case, approximately 40% of the total energy is ex- 
changed between the internal and center-of-mass degrees 
of freedom. The peaks of the center-of-mass kinetic en- 
ergy occur near quarter-periods of the harmonic motion. 
This surprising behaviour is due to the final coupling 
term in the ensemble Hamiltonian Eq. (|30|) . Alterna- 
tively, it may be attributed to the final nonlinear term in 
the modified Schrodinger equation l[7]). The nonlinearity 
explains why the dynamics are non-periodic. 

The distinctive hybrid dynamics, in which the center- 
of-mass kinetic energy is time-dependent, are the central 
result of this paper, and might be used to experimen- 
tally test the theory. To do so we must first find classical 
particles to experiment on. The predictions of the hy- 
brid theory give clues as to how such particles might be 
identified. 

There are no free parameters in the hybrid theory. 
Once the experimental parameters, such as masses, are 
known, the hybrid dynamics are quantitatively and ab- 
solutely determined. If those dynamics are not observed, 
the theory is refuted. 

Experimental signatures of hybrid dynamics can be 
found in position or momentum measurements. The 
most dramatic examples are perhaps in the dynamics of 
the center-of-mass kinetic energy En, shown in Fig. Ufor 
equal masses. Ensembles of momentum measurements 
allow En to be inferred, either directly if the center-of- 
mass momentum is measured, or indirectly if the indi- 
vidual particle momenta are measured. Momenta might 
be measured using scattering experiments. Both quan- 
tum and classical systems have constant center-of-mass 
kinetic energy, while only in hybrid systems does it vary 
in time. Hence the observation of a time- varying center- 



1.2 
1.0 

0.8 
0.6 
0.4 
0.2 
0.0 



E R /fiD , EjlfwjQ 



(a) 







Er 



12 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0. 



~ E R /haj , E,/haj 






























' E, 








; (b) 


: J 


I J 


\E R 


J 













10 



12 



FIG. 5: Hybrid oscillator center-of-mass kinetic energy Er 
(solid line) and internal energy Ei = E r + (V) (dashed line) 
versus dimensionless time t — io; M . (a) Dominant classical 
mass tn, = 20m q — 21/j,. (b) Dominant quantum mass 
m q — 20m x = 21[i. The vertical dotted lines are at quarter- 
periods. The contribution of the classical motion to the in- 
ternal energy has not been included. Initial conditions and 
other parameters are the same as for Figs. [5] and [3] 



of-mass kinetic energy would be evidence for hybrid dy- 
namics. 

Hybrid systems could also be distinguished from quan- 
tum or classical systems by ensembles of position mea- 
surements. For example, Fig. [3] shows that in the equal 
mass case, only in the hybrid dynamics does an initially 
zero covariance between the center-of-mass and relative 
coordinates evolve to differ from zero. 

Our numerical solutions may be related to specific 
experimental parameters through the reduced mass /j, 
and oscillator frequency uj^. The length scale x$ = 
(H/ iioj^) 1 / 2 and energy scale fuv^ follow. 

We next consider cases in which the classical and quan- 
tum particles have different masses. Fig. [5] shows the 
center-of-mass and internal energies in two cases, with ei- 
ther the classical, Fig. [5ja), or quantum, Fig.[5jb), mass 
dominant. Due to decoherence it is likely to prove eas- 
ier to keep a lighter object quantum-mechanical; hence 
the case with dominant classical mass, Fig. (5][a), is more 
likely to be the first to be realised experimentally. In 
that case, the exchange of energy between the center-of- 
mass kinetic energy En and the internal energy Ej is less 
than in the equal-mass case. The dynamics then better 
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FIG. 6: Hybrid oscillator center-of-mass kinetic energy En 
(solid line), internal energy Ei — E r + (V) (dashed line), and 
potential energy (V) (dotted line) versus dimensionless time 
t — tujp. Masses are equal. The total energy is E = 12.51fttJ M . 
The vertical dotted lines are at quarter-periods. The contri- 
bution of the classical motion to the internal energy has not 
been included. Initial conditions and other parameters are 
the same as for Figs. and 2] except K xx — K qq — 100/xq- 



approximate those of conventional quantum or classical 
physics, and the hybrid case is less distinctive. 

In the case with quantum mass dominant, Fig. [5jb), 
the exchange of energy exceeds that in the equal-mass 
case, and the signature of hybrid dynamics is clearer. 
As the ratio of the classical mass to the quantum mass 
becomes smaller, the peaks in the center-of-mass kinetic 
energy seen in Fig. EJb) become narrower. 

In the specific cases examined so far, the total en- 
ergy has been slightly greater than hw^. Fig. [5] shows 



the center-of-mass kinetic, internal and potential energies 
when the total energy is about 12^10^. This was achieved 
by reducing the initial coordinate variances, and hence 
increasing the initial momentum variances, by factors of 
100. When the potential energy is near zero, the center- 
of-mass kinetic energy spikes while the internal energy 
dips to conserve the total energy. Similar spiking be- 
haviour is also seen in Figs. U and [5] The spiking also 
occurs in the kinetic energy of the individual quantum 
particle, at the expense of a dip in the classical particle's 
kinetic energy. However, the total kinetic energy oscil- 
lates smoothly with no spiking. This behaviour is specific 
to the hybrid physics, and cannot be explained in terms 
of classical or quantum dynamics. 

In summary, we have applied the hybrid theory to the 
illustrative case of two harmonically coupled particles. In 
this case, we were able to use ansatzes to reduce the the- 
ory's partial differential equations to ordinary differential 
equations that are easy to solve numerically. The hybrid 
theory predicts different dynamics to classical physics or 
quantum physics. These differences are due to the dis- 
tinctive coupling of the relative and center-of-mass de- 
grees of freedom, and to the nonlinearity of the hybrid 
theory. Further work is required to develop an under- 
standing of the new hybrid dynamics. 

Applications of the theory to more realistic cases, such 
as nano-mechanical systems, will require numerical solu- 
tions of the theory's dynamical partial differential equa- 
tions, Eqs. (|5|6j) . An attractive feature of the hybrid 
theory is its lack of free parameters, which allows exper- 
iments to absolutely confirm or refute it. 
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